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ABSTRACT: 

Simultaneous  measurements  of  ocean  microstructure  and  sound  phase  shift 
from  a  stable  platform  in  BASS  STRAIT,  AUSTRALIA,  have  provided  new  relations 
between  the  statistics  of  the  medium  and  the  statistics  of  the  local  sound 
phase  speed  near  the  sea  surface  in  the  open  ocean.   Because  of  dispersion  due 
to  ambient  bubbles,  average  phase  speeds  in  the  frequency  range  15  to  100  kHz 
differ  as  much  as  2.5m/sec  from  the  accepted  3MHz  "precision"  velocimeter 
values  down  to  depths  of  7.6m  in  the  presence  of  wind  speeds  of  2  5-30  knots. 
These  differential  speeds  imply  average  bubble  volume  fractions  of  the  order 
of  10"'  with  standard  deviations  approximately  one-fifth  of  the  mean  value. 
The  differential  sound  speed  is  now  shown  to  increase  approximately  propor- 
tional to  the  wind  speed.   The  third  power  decrease  of  differential  speed  with 
increasing  depth  is  roughly  verified.   Under  these  experimental  conditions 
the  predominant  cause  of  the  local  phase  fluctuations  at  24.4  and  95.6  kHz  is 
shown  to  be  bubble  activity  rather  than  temperature  fluctuations.   At  24.4  kHz 
the  activity  is  the  random  change  of  number  of  bubbles.   At  a  frequency  such 
as  95.6  kHz,  where  there  is  a  large  resonant  bubble  population,  the  predomi- 
nant part  of  the  frequency  spectrum  of  the  sound  phase  modulation  is  shown  to 
be  caused  by  changing  bubble  radius  due  to  the  fluctuating  ocean  surface  wave 
height.   The  sound  phase  spectrum  mimics  the  wind  wave  spectrum  given  by 
Pierson  and  Moskovitz  to  two  octaves  beyond  the  frequency  of  the  peak  energy, 
at  which  point  the  surface  pressure  effect  has  dropped  low  enough  for  temper- 
ature fluctuations  to  take  over.   A  theory  is  presented  for  prediction  of 
these  microstructural  sound  phase  fluctuations  from  a  knowledge  of  the  surface 
wave  height  spectrum. 

This  task  was  supported  by  Code  302,  Naval  Ship  Systems  Command,  and  the 
Royal  Australian  Naval  Research  Laboratory. 

z:   //"  


TABLE   OF   CONTENTS 

Page 

I.  INTRODUCTION   --- 1 

II.  BACKGROUND -- - 1 

III.  THEORY 4 

A.  PHASE  SPEED  IN  A  BUBBLY  MEDIUM  -- 4 

B.  STATISTICS  OF  PHASE  SPEED  IN  A  FLUCTUATING  MEDIUM 10 

C.  THEORETICAL  BASIS  FOR  THE  EXPERIMENTAL  MEASUREMENTS 16 

IV.  EXPERIMENTAL  DETAILS 17 

V.  DATA  PROCESSING  - -- -23 

VI.  ANALYSIS  OF  RESULTS  - —24 

A.  AVERAGE  SPEED  OF  SOUND  AS  A  FUNCTION  OF  DEPTH —24 

B.  AVERAGE  SPEED  OF  SOUND  AS  A  FUNCTION  OF  SOUND  FREQUENCY— 28 

C.  AVERAGE  BUBBLE  FRACTION 30 

D.  FLUCTUATIONS  OF  BUBBLE  FRACTION  AND  SOUND  SPEED -30 

E.  TEMPORAL  CORRELATIONS  - — —36 

F.  SPECTRA —-45 

VII.  ACKNOWLEDGEMENTS —58 

VIII.  APPENDICES 

1.  DATA  PLANS  FOR  BASS  I,  II - - - -59 

2.  BATHY  THERMOGRAMS -—61 

3.  MATERIAL  VELOCITIES  CALCULATED  FROM  ELECTROMAGNETIC 

FLOWMETER  OUTPUT  DURING  BASS  II - -64 

IX.  REFERENCES 65 

X.  DISTRIBUTION  LIST  -- 67 


I.   INTRODUCTION 

The  objective  of  this  research  was  to  study  the  relations  between 
acoustic  phase  shifts  and  the  fluctuating  parameters  of  the  near  sur- 
face ocean.   Time  series  data  of  the  acoustic  phase  were  used  to  cal- 
culate the  dependence  of  the  standard  deviation  and  the  mean  value  of 
the  sound  speed  on  frequency,  range  and  depth.   Simultaneous  measurements 
were  made  of  the  fluctuating  temperature,  two  components  of  particle 
velocity  and  standard  acoustic  velocimeter  output. 

The  research  was  conducted  using  the  BHP/ESSO  Ocean  Drilling  Plat- 
form "Kingfish  A",  Lat.  38°,  39. 51  S,  Long.  148°,  08.7'  E,  Bass  Strait, 
southeast  of  Victoria,  Australia,  as  a  test  platform.   The  platform  is 
about  50  miles  from  shore  and  stands  in  water  approximately  250  feet 
deep.   Data  were  collected  during  two  visits  to  the  tower:   BASS  I 
(20-25  November  1972)  and  BASS  II  (12-18  February  1973). 

This  report  is  being  issued  jointly  with  the  Royal  Australian  Naval 
Research  Laboratory,  Garden  Island,  2000,  Australia. 

II.   BACKGROUND 

Following  the  early  work  of  Minnaert  (1933)  both  Germany  and  the 

2 
United  States  and  its  allies  conducted  extensive  theoretical  and  applied 

research  during  and  following  World  War  II  to  evaluate  the  importance  of 

ocean  bubbles  in  the  propagation  of  underwater  sound.   By  the  end  of  the 

1950s,  the  physical  understanding  of  bubble  action  in  the  presence  of  a 

3 
sound  wave  was  secure  and  in  the  1960s  research  spread  to  many  adjacent 

areas  including  chemical  engineering,  nuclear  physics,  and  oceanography 

and  even  meteorology  . 

1. 


The  question  of  whether  there  are  bubbles  in  the  ocean,  without  the 
intrusion  of  man-made  vehicles,  has  remained  controversial.   Those  who  tend 
to  doubt  the  existence  of  ambient  bubbles  in  the  sea  quote  the  laboratory 
experiments  of  Liebermann  ,  who  demonstrated  the  rapid  loss  of  bubbles  due 
to  buoyancy  and  diffusion.   On  the  other  side  of  the  controversy,  there  have 
been  measurements   showing  that  bubbles  are  produced  by  wind  waves  and  decay- 
ing matter  on  the  ocean  floor  .   Marine  biologists  have  no  doubts  about  the 
continual  generation  of  bubbles  in  the  sea  by  photosynthesis,  and  marine 
life.   The  root  question,  of  course,  is  whether  those  bubbles  that  are  pro- 
duced by  physical  and  biological  activities  in  the  sea  are  so  numerous  that 
in-situ  measurements  will  always  detect  ambient  bubbles  and,  more  directly 
for  the  underwater  acoustician,  the  question  is  whether  ambient  bubbles  are 
numerous  enough  to  influence  sound  propagation. 

The  Ocean  Physics  Laboratory  at  the  Naval  Postgraduate  School  conducted 

a  student  research  program  to  determine  the  excess  attenuation  and  scatter 
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from  near  surface  bubbles  from  1963  through  1965.   This  early  work  *  was 

summarized  in  a  recent  publication   which  showed  the  results  of  photograph- 
ing bubbles  and  of  in-situ  acoustical  attenuation  and  scattering  measurements 
at  the  Naval  Underwater  Center's  oceanographic  tower  in  San  Diego,  California, 
and  in  Monterey  Bay.  Measurements  were  made  from  the  surface  to  depths  of 
45  ft.   Bubble  volume  fractions  of  the  order  of  10    or  greater  in  a  one 
micron  radius  increment  were  typical.   The  dominant  sound  frequencies  that 
would  be  affected  were  in  the  measurement  range,  above  about  25  kHz,  but  theory 
indicates  that  significant  dispersion  of  sound  speed  should  also  be  found  for 
frequencies  lower  than  the  measured  ones. 

The  Ocean  Physics  Laboratory  set  out  to  measure,  directly,  the  change  of 


sound  speed  with  sound  frequency,  implied  by  the  previous  determinations  of 
attenuation  and  scatter  near  the  sea  surface.   Results  of  measurements   from 
the  NUC  tower  at  San  Diego  in  October  1971,  showed  unmistakable  evidence  of 
dispersion  of  the  speed  of  sound.   Deviations  of  the  propagation  speed  with 
frequency  were  of  the  order  of  meters/second  greater  or  less  than  the  "theo- 
retical" speed  calculated  from  the  traditional  parameters  of  temperature, 
salinity  and  depth.   The  standard  deviations  of  the  speed  measurements  were 
much  less  than  the  deviations  of  the  measured  speed  from  the  theoretical 
speed.   Bubble  populations  were  identifiable  at  depths  4.3m  and  9.3m,  which, 
when  converted  to  surface  values  on  the  basis  of  constant  mass  and  isothermal 
transport,  represent  62  micron  (52  kHz)  bubbles  and  100  micron  (32  kHz) 
bubbles.   The  power  spectral  density  of  the  time-varying  sound  speed  peaked 
at  frequencies  less  than  0.5  Hz,  suggesting  that  bubbles  were  entrained  at 
orbital  frequencies  associated  with  the  ocean  surface  wave  system. 

During  1972,  our  group  mounted  three  simultaneous  attacks  in  order  to 

understand  better  the  dynamics  of  bubble  populations  in  the  upper  ocean.  Trans- 

12 
port  analysis   of  the  effects  of  buoyancy,  gas  diffusion  and  vertical  velo- 
cities of  the  medium  showed  that  bubbles  which  originate  at  the  surface  would 
require  large  downward  convective  velocities  to  persist  at  depth;  large 
bubbles  which  originate  in  the  sediment  rise  to  the  surface  with  little  evi- 
dence of  gas  diffusion;  an  exponentially  decaying  volume  source  (maximum  at 
the  surface)  would  show  ambient  distributions  that  mimic  the  source  function, 
if  the  bubbles  are  large,  and  would  show  much  greater  variation  from  the 

source  distribution  as  a  result  of  gas  diffusion,  if  the  bubbles  are  small. 

13 

A  study   of  the  statistics  of  fluctuations  of  sound  speed  due  to  bubble 

populations  showed  that   the  maximum  fluctuations  due  to  changes  in  number  of 


3. 


bubbles  would  occur  at  frequencies  displaced  from  the  resonant  frequency  by 


amoun 


ts  fD/(l  ±  S)2.   The  standard  deviation  would  be  greater  at  the  lower 
R 


frequency  than  at  the  upper  one.  fR  =  resonant  frequency,  6  =  damping  constant 

14 
An  experimental  study   of  sound  speed  at  the  Naval  Underwater  Center 

oceanographic  tower  in  June  1972,  covered  the  extended  frequency  range  15  to 

150  kHz  in  increments  of  approximately  3.5  kHz  at  a  depth  of  3.1  ±  0.3  meters. 

Both  average  speed  and  standard  deviation  of  speed  were  determined  at  each 
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frequency.   An  analysis   of  the  experimental  data  showed  that  the  statisti- 
cal theory  predictions  of  frequencies  of  maximum  fluctuation  of  sound  speed 
agreed  with  the  experimental  measurements  for  at  least  four  identifiable 
bubble  populations  of  surface  equivalent  radii  54,  76,  105  and  200  microns, 
resonant  at  surface  equivalent  frequencies  of  60,  43,  31,  16  kHz,  respec- 
tively.  Correlations  between  sound  phase  and  ocean  surface  height  revealed 
that  there  are  strong  fluctuations  in  the  sound  speed,  in  phase  with  the 
surface  height  fluctuations,  at  sound  frequencies  that  correspond  to  the 
resonant  frequencies  of  the  peak  bubble  concentrations. 

The  purpose  of  the  BASS  experiments  was  to  repeat,  in  the  waters  of  the 
South  Pacific,  some  of  the  San  Diego  experiments  and  to  look  for  the  statis- 
tics of  the  correlations  with  speed  fluctuations  and  its  temperature  micro- 
structure  and  medium  particle  velocity  as  a  function  of  depth  and  sound 
frequency. 

III.   THEORY 

A.   Phase  Speed  in  a  Bubbly  Medium 
1.   Bubbles  of  One  Size  Only 

The  presence  of  bubbles  in  a  water  medium  affects  the  speed  of  sound 
because  of  the  changed  compressibility.   Although  there  is  also  a  potential 

4. 


effect  on  the  speed  caused  by  the  reduced  density  of  the  medium,  we  are 
interested  only  in  the  ocean  medium,  for  which  the  density  change  can  readily 
be  shown  to  be  negligible  for  realistic  environments. 

The  dependence  of  the  speed  on  the  compressibility  and  density  is  given 
by 

2  -   1  /IN 

c   =7F  (1> 

The  compressibility  of  the  bubbly  medium,  K,  is  made  up  of  a  part  due  to  the 
water,  KQ ,  and  the  part  involving  the  change  in  volume  of  the  bubbles,  Kx , 

K  -  K  +  K 
o    i 

The  compressibility  of  the  bubble-free  water  is  most  simply  expressed  in  terms 
of  the  usual  value  of  the  speed  of  sound  through  water,  c  ,  and  the  density, 

<V 

The  speed  through  the  bubble-free  water  can  be  obtained  from  a  simplified  form 
of  Wilson's  equation 

c0  =  1449.  +  4.6t  -  0.055t2  +  0.017d  (2a) 

where  t  =  temperature  [°C],  d  -  depth  [meters],  and  salinity  is  assumed  35  ppt. 

The  more  sensitive  compressibility,  due  to  the  bubbles,  is  given  by 

2 
AV   H&y   NS  |  NS 


P     P     P     muu  ( -1  +  _£)  +  J  _ 

LN      uu27     U)m  J 

where     N  =  number  of  bubbles  of  radius,  a 

Av  =  change  in  volume  for  each  bubble  in  incident  field,  P 

2 
S  ■  4TTa  =  surface  area  of  each  bubble 

and     I  =  complex  magnitude  of  radial  displacement  of  bubble  surface. 
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The  compressibility,  K,  is  written  in  complex  notation  because  of  its  depen- 
dence  on  the  radial  displacement,  ? >  ,  which  is  itself  complex.   Simplifica- 
tion  occurs  if  we  now  insert  §_  and  also  define 

Z  =  fR/f  -   %/uu  (4) 


and 


Then, 


6  =  R/uu 
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N4TT3  [Z   ~  1  ~  j&] 


2  r/^2 

p  w  L(z 


-  1)' 


62] 


(5) 


(6) 


and 


K  =  K  +  K  = 


'oc, 


4naNLz  -  1  ~j6j 
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The  expression  for  the  speed  of  sound  in  the  bubbly  medium  is  now 
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and  the  real  part  of  c  is 
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Before  interpreting  our  result,  it  is  useful  to  put  the  bubble  measurement 
in  terms  of  the  fraction  of  gas  in  bubble  form,  U(a),  defined  for  bubbles 
of  radius,  a,  by 

U(a)  =  [N(a)][jTTa3] 

Then, 
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It  is  of  interest  to  look  at  the  special  cases  for  frequencies,  (a)  much 
less  than  the  resonant  frequency  of  the  bubbles,  (b)  much  greater  than  the 
resonant  frequency,  and  (c)  close  to  the  resonant  frequency.   For  simplic- 
ity, the  reduction  is  carried  out  for  standard,  sea  level  conditions. 


(a)   f  «  fR,   Z  »  1 


[-*] 


c  =  c0|i -rr  (ii) 


-2 

then,  since  at  sea  level  al<_  =  1.36  x  10   ,  the  speed  reduces  to 

c  =  c  [l  —  8.1  x  10  U(a)]  and  the  value  is  found  to  be  less  than  c_ ,  inde- 
pendent of  the  particular  frequency,  f,  and  dependent  only  on  the  fractional 
volume  of  gas  to  water,  U(a). 

(b)    f  »  fR,   Z  «  1 

At  sea  level,  the  speed  has  the  value 


1  +  8.1  x  103  U(a)  — - — =• 
L  1  +  62J 


c  =  co|  1  +  8.1  x  10  U(a)      0  j  (12) 

which  is  greater  than  the  bubble-free  value  c  ,  dependent  on  the  particular 
frequency  and  damping  constant,  reaching  the  value  c  only  as  the  frequency 
approaches  infinity. 


(c)   f  «  fR,  Z  «  1 

The  speed  in  the  bubbly  water  is  c  at  the  resonant  frequency. 
It  rises  rapidly  to  a  peak  at  frequencies  slightly  greater  than  fR,  and  drops 
rapidly  to  a  minimum  at  frequencies  slightly  less  than  f  . 

K 

The  detailed  behavior  of  the  speed,  near  the  resonant  frequency,  is 


7. 


shown  in  Fig.  1,  which  has  been  drawn  for  two  cases  of  different  numbers  of 
bubbles.   The  change  in  speed  has  been  presented  in  terms  of  the  shift  in 
speed  6c  =  (c  —  c  )  at  frequency  ratio  f/fR'   It  is  interesting  to  observe 
that  the  maxima  and  minima  of  the  resonant  swings  occur  at  frequencies  given 
approximately  by  f  =  fR(l  ±  6R/2).   That  is,  for  6R  =  0.10,  the  peak  and 
trough  of  the  resonance  are  6R/2  ~  57«  above  and  below  the  resonant  frequency. 
The  magnitudes  of  the  speed  at  these  points  are  above  and  below  the  speed 
in  bubble-free  water  by  an  amount  approximately  rr-  times  the  value  at  the 
zero  frequency  limit  given  by  equation  (11). 

Study  of  equation  (10),  which  led  us  to  Figure  1,  shows  that  the  speed 
of  sound  is  significantly  altered  even  at  frequencies  far  from  the  reson- 
ance frequency.   For  example,  for  the  case  of  6_  -   0.1,  the  speed  at 

3 
f/f  =  2  is  still  approximately  0.33  (8.1  x  10  Uc  )  above  the  bubble-free 
R  o 

value  which  exists  as  f  approaches  infinity.   And  at  f/fD  =  0.5,  the  lowered 

K 

3 
value  of  the  speed  is  0.33  (8.1  x  10  Uc  )  less  than  the  asymptotic  value 

achieved  as  the  frequency  approaches  zero. 

2.   Bubbles  of  many  Sizes. 

The  generalization  of  our  results  to  a  cluster  of  bubbles  of 
various  radii  is  accomplished  by  including  additional  terms  in  equation  (7). 
It  is  easy  to  handle  the  situation  for  a  broad  range  of  radii,  such  as 
exists  in  the  ocean,  by  replacing  N(a)  by  n(a)da  =  number  of  bubbles  per 
unit  volume  with  radius  between  a  and  a  +  da.   Similarly,  U(a)  is  replaced 
by  u(a)da  to  represent  the  volume  fraction  of  gas  in  bubble  form  to  liquid. 
Then,  instead  of  K  ,  as  in  equation  (7),  we  write 

41raln.(al)dal   [z^-  1-  j61R] 
^"A      V2         [(V-l)2^2]  (13) 
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Because  all  contributions  to  the  compressibility  are  very  small  quantities, 
they  add  linearly  and  the  speed  of  sound  in  the  bubble  region  can  be  written 


Re{c}  =  c 

Au/         O 


1  - 


u(a,)Z.2(z.2  -  lW 
i  1  \  1     /   l 


3  f      "<ai>zAzi  -1K 


(14) 


The  effect  of  a  mixture  of  bubble  sizes  around  a  peak  population  at  reson- 
ance frequency,  f  ,  is  to  smear  the  dispersion  curve  for  the  speed  so  that, 
although  the  magnitudes  of  the  deviation  from  the  bubble  free  value  are 
increased,  the  frequency  range  between  the  peak  and  the  trough  is  also  in- 
creased beyond  the  value  that  we  calculated  in  Section  A  for  the  single 
bubble  size. 

B.   Statistics  of  Phase  Speed  in  a  Fluctuating  Medium 

The  near-surface  region  of  the  ocean  is  a  medium  that  shows  spatial  and 
temporal  fluctuations  of  the  temperature,  salinity,  pressure,  and  three  com- 
ponents of  particle  velocity.   Sound  propagation  takes  place  in  such  a 
medium  only  at  the  cost  of  sound  phase  and  amplitude  fluctuations  along  the 
acoustic  path.   The  phenomena  have  been  studied  theoretically  under  rather 
major  simplifying  assumptions   .   Recent  experimental  research  has  extended 
these  studies  to  the  near-surface  region  of  the  ocean   ,  in  order  to  search 
for  the  real  world  correlations  between  the  sound  fluctuations  and  the 
fluctuations  of  the  medium. 

Fluctuations  in  the  speed  of  sound,  measured   near  the  ocean  surface, 
have  been  shown  to  have  a  Gaussian  probability  density  function  except  in 
the  regions  of  large  differential  speed  of  sound,  where  the  PDF  was  double- 
peaked.   The  temporal  correlation  time  for  phase  was  about  the  same  as  for 
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the  oceanographic  variables.  The  power  spectral  density  of  the  time-varying 
phase  showed  evidence  of  bubbles  entrained  at  orbital  frequencies  associated 
with  the  surface  wave  system. 

Work  at  the  Naval  Postgraduate  School   has  also  revealed  that  orbital 
motion  due  to  surface  wave  action  shifts  inhomogeneities  of  temperature, 
salinity,  and  bubbles,  and  causes  phase  changes  of  sound  propagated  near  the 
surface.   Phase  modulation  envelopes  were  found  to  have  smaller  correlation 
times  and  spatial  correlation  lengths  in  the  mixed  layer  than  in  the  thermo- 
cline.   Correlation  times  were  greater  for  propagation  in  the  vertical  than  in 
the  horizontal  directions  in  the  thermocline.   Anisotropy  in  the  thermocline 
was  also  indicated  by  the  variance  of  phase  fluctuation  being  greater  for 

horizontal  than  vertical  propagation. 

14 
Cross-correlations  between  ocean  surface  height  and  sound  phase   have 

been  shown  to  be  very  large  at  the  sound  frequencies  where  dispersion  of  the 

speed  indicated  large  populations  of  bubbles. 

For  the  purpose  of  this  report  it  will  not  be  necessary  to  look  at  all 

of  the  details  of  the  mechanisms  for  sound  speed  fluctuations.   The  problem 

can  be  discussed  adequately  if  we  simply  point  out  that  the  variance  of  the 

speed  fluctuation  can,  in  principle,  be  due  to  any  of  the  medium  descriptors 

mentioned  in  the  first  paragraph  of  this  section.   However,  realistically, 

fluctuations  in  salinity  and  density  can  readily  be  shown  (from  eq  (2a)  to 

be  negligible.   Furthermore,  it  is  possible  to  subtract  out  the  effects  of 

motion  of  the  medium  either  by  a  double  path,  as  in  the  case  of  commercial 

velocimeters,  or  by  subtracting  out  the  effect  of  the  motion  indicated  by  a 

velocity  flowmeter,  as  we  have  done.   It  is  the  latter  technique  that  has 

been  used  in  this  experiment.  This  leaves  fluctuating  bubbles  and  temperature 
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as  the  major  sources  of  sound  speed  variation. 

It  is  possible  to  calculate  the  variance  of  the  speed  fluctuations  due 
to  temperature, simply  and  effectively  enough  for  the  purposes  of  this  study, 
by  linearizing  equation  (2a)  and  rewriting  it  as 

c0  -  cref  +  bt  (15) 

For  water  of  temperature  19°  C,  the  linearization  parameters  are: 

c  r,  -    1471.25  m/sec 
reD  =    2.66  fc)"1 

The  variable  effect  of  bubbles  is  due  to  (a)  changes  of  total  gas  in 
bubble  form,  U,  or  (b)  changes  of  the  bubble  distribution  with  radius,  or 
(c)  perturbation  of  the  resonant  frequency  of  a  predominant  population  at 

the  frequency  of  the  study.   The  source  of  fluctuation  (a)  has  been  con- 

13 
sidered  carefully  '  and  it  has  been  shown  that  the  maxima  occur  at  frequencies 

given  by  f  =  fR(l  ±  ^r/2)  for  resonant  frequency,  f  ,  with  damping  constant, 

V 

There  is  little  information  about  preferential  sizes  of  bubbles  for 
specific  bubble-making  processes.   However,  there  is  no  reason  to  assume  that 
the  average  distribution  of  bubble  numbers  with  radius  at  a  given  depth  and 
time  of  day  will  change,  except  by  rather  small  amounts  during  a  propagation 
experiment.  We  will  assume  that  fluctuation  source  (b)  .is  negligible. 

Finally,  (c)  perturbation  of  the  resonant  frequency  can  have  a  profound 
effect  on  the  variation  of  the  speed  of  sound  and  it  appears  that  the  change 
of  water  pressure  (due  to  surface  wave  action)  at  the  level  of  the  experiment 
can  modulate  the  resonant  frequency  as  given  by  the  equation 


■  -—   f 

"R   2na  V 


3YP 


P 
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where    P  =  ambient  pressure 
o 

p   =  ambient  density 

o  J 

and  thereby  cause  the  speed  of  sound  to  move  up  and  down  the  steep  part  of 
the  speed  curve,  as  shown  in  Figure  1.   If  this  action  occurs  during  propa- 
gation of  a  frequency  that  corresponds  to  a  major  bubble  population,  a  very 
large  variance  of  sound  phase  can  occur. 

The  effect  of  bubbles  in  changing  the  speed  of  sound  can  also  be  lin- 
earized if  we  restrict  our  interest  to  one  frequency  at  a  time.   Then,  in  a 
bubbly  medium 

c  =  c0  [1  +  6c(U,f)]  (17) 

where    6c(U,f)  is  the  increment  in  speed  due  to  bubbles.   Combining  equa- 
tions (15)  and  (17),  since  both  bT  and  6c(U,f)  are  very  small  compared  to 
unity,  the  expression  for  the  speed  of  sound  can  be  written 

c  ~  c   ,  +  bT  +  6c  (18) 

ref 

from  which  the  variance  of  the  speed  is 

Var  {c}  =  b2  Var  [t]  +  Var  {6c}  (19) 

Experimental  values  of  the  variance  of  6c  can  be  used  in  two  ways: 

(a)  At  very  low  frequencies,  equation  (11)  leads  to 

9c2 
Var   {6c}  =        ,  %      Var   {u}         ,        f  «  fR    ,  (20) 

4a2kR2  R 

so  that  the  variance  of  changes  of  total  fraction  of  bubble  volume,  Var  lUj, 
can  be  calculated. 

(b)  At  frequencies  that  show  1007o  cross-correlation  of  sound 
phase  with  the  surface  wave  height,  it  can  be  assumed  that  the  high 
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correlation  is  completely  due  to  resonant  frequency  fluctuation  caused  by 
pressure  and  radius  perturbations  defined  by  equation  (16)  near  Z  =  fD/f  =  1 
A  change  of  fR  at  constant  f  can  be  understood  by  considering  the  speed 
change  shown  in  Figure  1.   At  f/fR  =  1»  there  is  a  maximum  rate  of  change 
of  c  with  change  of  the  ratio.   For  bubbles  beyond  the  peak  and  trough  of 
the  dispersion  curve,  the  effect  of  change  of  t/f.      is  very  much  less,  and 
goes  in  the  opposite  direction.   A  positive  100%  cross-correlation  with 
wave  height  would  suggest  that  phase  changes  near  Z  =  1  override  changes  in 
the  opposite  direction  in  the  outer  regions  of  the  dispersion  curves  of  the 
much  large  and  much  smaller  bubbles. 

We  consider  the  region  near  Z  =  1  by  defining 

Z  -   f_/f  =    1  +  A  where  A  «   1  (21) 

K 


V  3YP  /  P 

A  =  ~nrtt —  l    for  z « i  (22) 

For  the  slow  changes  accompanying  the  ocean  wave  system  we  can  assume  the 
isothermal  relation 

da/dP=-^  . 
Differentiation  of  the  above  equation  then  gives  the  rate  of  change 

f  ■  £  (1  +  A>  <23> 

Integration  between  limits 


\  10m/ 


Ax  =  Zl  -   1  =  0  at  P1  =  P  I  1  + 
*2=*  atr2-Ps(1  +  T^  +f^ 
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yields 


*-* 


ffc 


[.    h,m/10m 
10m_ 


(24) 


where  P  ■  surface  ambient  pressure  ■  1  atmosphere 
d,m  =  experiment  depth  in  meters 
h,m  =  wave  height  in  meters 

The  expansion  for  small  increments  is  used 


&i(l  +  x)  =  x  -  j   x2  + 


.         h,m/10m 
where  x  =  — * - — 

10m 


to  give 


A  w| 

6 


[h,m/10m 
1  +  ^H 
10m  J 


(25) 


Using  the  approximation  Z  =  1  +  A  and  taking  the  variance  of  both  sides  of 
equation  (10),  we  get 


Var  {6c (A)}  =  Var  i  c 


1  - 


to  order  A  «  1.    (26) 


YPfl  (4A 


2     2 
If  we  assume  that  4A  «  6   (e.g.,  much  less  than  a  5%  variation  in  f^/f 

R 


for  a  60  micron  bubble)  and  use  equation  (24),  we  get 

r    r 


Var   {6c (A)}  =  Var  , 


or 


"■fcW  = 


0 

V        «- 


0.83   Uc 


1  - 


0.83 


h,m 
10m 


.     2    2c2     ...         d  ,m 
kRa   6      (1+T0m" 


(vK 


1  +  TS5>2 


h/10 


(27) 


(28) 


where  O""..  (A)  =  rms  value  of  fluctuations  in  sound  speed  due  to 

perturbations  of  bubble  resonant  frequency. 
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°   L/ir,  =  tne  dimensionless  fraction,  rms  height  of  surface  wave 
h/  iU 

in  meters/10  meters. 
Since  we  are  dealing  here  with  very  slowly  varying  surface  heights 
(F  <  1.  Hz)  the  change  of  resonant  frequency  and  radius  of  the  bubble  will 
occur  with  negligible  phase  lag.   Therefore,  equation  (25)  can  be  viewed 
from  the  point  of  view  of  a  sinusoidal  change  of  height,  h,m  =  H,m  sin  Qt 
where  Q  =  2ttF 

and  F  =  a  component  of  the  surface  wave  spectrum  causing  a 
simultaneous  change  of  A  in  equation  (25)  and  of  the  differential  sound  speed 
without  phase  change  at  the  same  driving  frequency.   In  the  case  of  a  spec- 
trum of  surface  height,  equation  (28)  can  be  rewritten 


Sc(W0-gpS  (28a) 


10my 


C.   Theoretical  Basis  for  the  Experimental  Measurements 

The  physical  form  of  the  acoustic  sound  path  is  shown  in  Figures  3  and 
4.   Sound  from  the  piston-like  transducer  propagates  through  the  frame  and 
the  sound  pressure  is  measured  by  a  series  of  very  small  "identical"  hydro- 
phones along  the  path.   For  a  plane  wave,  or  in  the  far  field  of  the  piston, 
the  phase  shift  along  the  distance,  x,  between  two  hydrophones  is  measured 
by  counting  the  integral  number  of  wave  lengths,  n,  and  determining  the 
residual  phase  angle,  0,  which  represents  the  instantaneous  fractional  wave 
length,  0/360.   Then  the  instantaneous  speed  is  calculated  by  the  product  of 

frequency  by  wave  length,  using 

fx 
C  =  n  +  0/360  (29) 

The  assumptions  in  the  derivation,  as  well  as  the  accuracy  of  the  component 

measurements,  are  discussed  in  reference  11. 
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IV.   EXPERIMENTAL  DETAILS 

A  photograph  of  Kingfish  A  Tower  is  presented  in  Figure  2.   A  rectan- 
gular aluminum  pipe  frame  of  height  six  feet  and  length  eighteen  feet,  was 
used  to  mount  the  transducers,  hydrophones,  flowmeter,  velocimeter  and  ther- 
mister  (Figs.  3  and  4).   To  prevent  displacement  of  the  hydrophones,  the 
hydrophone  mounting  wires  were  put  under  tension  and  isolated  from  the  frame 
by  rubber  shock  mounts.   To  minimize  any  reflections  from  the  frame,  selected 
segments  were  covered  with  acoustic  absorbent  rubber  (SOAB).   See  also  ref- 
erences 11  and  14. 

The  rectangular  frame  was  hung  so  that  the  sound  source  and  receivers 
were  in  a  horizontal  line.   The  method  of  attachment  to  a  nominal  45   steel 
cable  is  shown  in  Figure  4.   The  steel  cable  had  one  end  held  to  the  ocean 
floor  by  an  anchor,  the  other  end  was  fixed  to  the  tower.   The  steel  cable 
and  polypropylene  attachment  lines  can  be  seen  in  Figure  2  and  to  the  right 
in  Figure  3  as  they  enter  the  water. 

At  the  center  of  one  end  of  the  frame,  pointing  in  the  horizontal  (X) 
direction,  an  USNRD  F-33  piston  transducer  was  mounted.   On  the  acoustic 
axis  of  this  transducer,  four  Atlantic  Research  LC-10  hydrophones  were 
placed  at  source  separations  shown  in  Table  1. 

TABLE  1. 

SOURCE  TO  HYDROPHONE  SEPARATIONS  (METERS) 

BASS  I  BASS  II 

LC  10-1  1.047  1.065 

LC  10-2  0.998  0.988 

LC  10-3  2.002  1.990 

LC  10-4  1.355  0.986 
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The  transmitting  instrumentation  (Figure  5)  consisted  of  a  General  Radio 
coherent  decade  frequency  synthesizer,  type  1162-A,  generating  a  stable  sinus- 
oidal output  of  two  volts  rms,  selectable  in  frequency.   The  signal  was  ampli- 
fied by  a  Hewlett-Packard  Power  amplifier,  467A,  to  voltage  levels  consistent 
with  undistorted  signals  through  the  receiving  instrumentation,  and  then 
impressed  across  the  transducer. 

The  receiving  instrumentation  (Figure  5)  consisted  of  four  LC-10  hydro- 
phones with  flat  (-HOdB  ±  3  dB  re  1  volt/microbar)  receiving  responses  from 
10  kHz  to  150  kHz.   The  signals  from  the  hydrophones  were  amplified  30  dB 
by  NUS  pre-amplifiers  (model  2010-30)  before  passing  through  180  feet  of 
water-proof  shielded  cable.   Two  selected  signals  were  then  bandpass  filtered 
by  digitally  tuned  variable  Krohn-Hite  model  3322  filters  (down  24  dB/octave 
outside  of  band)  set  with  the  center  frequency  equal  to  the  transmitted 
frequency  ±  300  Hz.   The  filtered  signals  were  then  fed  into  a  Dranetz  model 
305  phase  meter.    The  phase  meter  measured  the  phase  angle  between  the  two 
ac  voltages  and  provided  a  dc  voltage  output  (10  ms  response  time)  propor- 
tional to  the  difference  of  the  phase  angles  at  10  mv  per  degree.   The  dc 
voltage  was  FM  recorded  at  3.75  ips  on  a  Precision  Instrument  model  2101, 
seven  track  magnetic  tape  recorder. 

The  thermister  was  mounted  on  the  frame  midway  between  the  first  two 
hydrophones.   The  dc  voltage  output  was  FM  recorded  at  0.5  v  per  degree  C. 

The  reference  sound  speed  was  obtained  by  an  NUS  model  2030  -  103  -  12 
velocimeter  operating  at  3MHz  transmit  frequency.   The  frequency  output  was 
time  averaged  and  manually  recorded.   The  frequency  output  was  also  demodu- 
lated to  produce  a  time  varying  dc  voltage  which  was  FM  recorded  at  20.48 
mv  per  m/sec  sound  speed. 
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Two  orthogonal  horizontal  components  of  the  water  particle  velocity 
were  determined  by  the  mounting  of  an  EPCO  water  current  meter  such  that  the 
velocities  measured  corresponded  to  the  direction  of  acoustic  propagation 
of  the  USNRD  F-33  transducer  and  the  transverse  or  Y  direction  (BASS  II). 
These  two  outputs  were  also  FM  recorded  at  1.640  volts  per  M/sec. 

V.   DATA  PROCESSING 

Data  processing  was  essentially  the  same  as  described  in  previous  Naval 
Postgraduate  School  reports  (Refs.  11,14).   See  Figure  5. 

A  series  of  ten  minute  and  four  minute  individual  records  makes  up  the 
analog  recorded  time  varying  sound  phase  difference,  water  current  components, 
temperature,  and  velocimeter  data.   The  approach  used  was  to  digitize  the 
analog  data  and  use  a  computer  to  carry  out  the  statistical  correlation  and 
spectral  analysis. 

The  recorded  data  were  digitized  using  the  Royal  Australian  Naval  Research 
Laboratory's  PDP-11  computer.   The  sampling  frequency  was  6.4  Hz,  giving  a 
Nyquist  frequency  of  3.2  Hz.   The  sampling  interval  was  0.15625  seconds  and 
the  frequency  resolution  was  0.05  Hz.   The  analog  signals  were  low-passed 
filtered  at  3.2  Hz,  scaled  and  written  on  magnetic  disk.   The  horizontal  beam 
(X)  direction  component  of  water  current  was  scaled  to  meters  per  second. 
For  each  instantaneous  value  of  phase,  the  instantaneous  value  for  the  speed 
of  sound  was  calculated,  in  meters  per  second,  and  instantaneous  value  for 
water  current  in  the  X  direction  was  algebraically  subtracted  to  correct  for 
apparent  sound  velocity  due  to  the  motion  of  the  medium  with  respect  to  the 
sound  field. 

The  mean  values,  variances  and  standard  deviations  were  then  determined. 
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The  mean  speed  was  then  compared  with  the  mean  speed  values  given  by 
the  NUS  Velocimeter. 

The  temporal  autocorrelations  were  computed  and  plotted  using  a  maximum 
time  lag  of  20  seconds.   The  temporal  cross  correlation  was  also  computed 
between  sound  speed  and  water  cross  current,  and  sound  speed  and  temperature. 
From  these,  the  auto-spectrum  was  computed. 

VI.   ANALYSIS  OF  RESULTS 

A.   Average  Speed  of  Sound  as  a  Function  of  Depth 

The  equipment  used  for  the  phase  comparison  determination  of  the  speed 
of  sound  is  capable  of  yielding  values  with  an  absolute  error  of  less  than 
one  meter  per  second   ,  provided  that  corrections  are  made  for  ocean  currents 
and  the  small  phase  difference  between  the  hydrophone-preamplifier  electrical 
responses.   Column  3  of  Table  2  presents  the  values  of  the  average  sound 
speed  (c)  determined  by  our  phase  comparison  method,  corrected  for  ocean  cur- 
rent, at  two  frequencies,  24.4  and  95.6  kHz  (BASS  II).   Hydrophone  separa- 
tion was  0.988  m.   Column  4  presents  the  concurrently  read  and  averaged 
values  of  the  mean  speed  (c0)  given  by  the  sound  velocimeter.   In  column  5, 
we  take  the  difference  between  column  3  and  4  in  order  to  arrive  at  the 
differential  velocity,  Ac  =  <c>  —  (c  > .   These  values  are  due  to  bubbles 
plus  a  correction  for  the  previously  mentioned  hydrophone  electrical  phase 
difference.   A  complete  frequency  run  at  depths  great  enough  to  validate  the 
assumption  of  a  no-bubble  region  at  13.7  m  was  not  concluded  because  of  the 
failure  of  the  tape  recorder.   However,  from  the  evidence  of  the  standard 
deviation  of  sound  speeds  shown  in  Table  3,  it  is  quite  convincing  that  the 
depth  of  13.7  m  can  be  used  to  represent  the  bubble-free  medium,  and, 
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therefore,  to  correct  for  the  hydrophone-preamplifier  electrical  phase  shift. 
This  has  been  done  to  arrive  at  column  6  of  Table  2,  the  corrected  differen- 
tial speed  of  sound  due  to  bubbles,  6c  -  Ac  —  Ac.  _  7.   Finally,  column  7 
gives  the  average  bubble  fraction,  U,  calculated  from  equation  (11),  using 
average  values. 

The  frequency  24.4  kHz  was  selected  because  it  was  the  lowest  frequency 
available  with  the  35  dB  S/N  ratio  required  for  accurate  phase  information. 
The  values  of  the  corrected  differential  speed  become  more  negative  (down  to 
6.1  M)  as  depth  is  increased,  presumably  because  there  is  a  peak  bubble  popu- 
lation just  above  24  kHz  at  surface  level  which  increases  its  resonant  fre- 
quency to  higher  values  with  increasing  depth.   A  strong  population  sweeping 
by  the  sound  frequency  of  24.4  kHz  would  cause  larger  negative  values  of  6c. 
Finally,  however,  by  7.6  m  depth,  the  decreasing  numbers  of  bubbles  becomes 
the  dominant  factor  and  the  speed  approaches  the  no-bubble  value. 

The  frequency  95.6  kHz  presents  a  continuously  decreasing  value  of  6c 
as  depth  increases.   Reference  14  shows  that  there  is  a  trough  of  sound  speed 
at  this  frequency,  caused  by  the  larger  number  of  bubbles  at  higher  resonant 
frequency.   As  the  depth  is  increased,  existing  large  bubbles  have  increas- 
ingly high  resonant  frequencies  which  produce  a  positive  increment  to  add  to 
the  speed  of  sound,  thereby  reducing  the  negative  corrected  fifferential 
speed  as  the  depth  is  increased.   The  smaller  number  of  bubbles  at  greater 
depths  acts,  also,  to  bring  the  value  of  6c  closer  to  zero. 

The  ratio  of  the  average  change  in  sound  speed  ascribable  to  bubbles, 
6c,  to  the  average  speed  in  non-bubbly  water  (c  )  (equivalent  to  the  veloci- 
meter  reading  at  3  megahertz)  is  plotted  in  Figure  6  as  a  function  of  the 
pressure  parameter,  p  =  1.33(1  +  0.1  z) .   In  reference  15  it  was  shown  that 
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FI6  6   Variation  of  the  relative  change  of  speed  due  to 
bubbles,  6c/ (c  y  as  a  function  of  the  pressure  parameter, 
p  =  1.33/(1+0.1  Z)  where  Z  =  depth  in  meters.   The  data  of  this 
experiment,  at  24.4  and  95.6  kHz,  measured  during  winds  speeds 
of  25-30  knots,  are  compared  with  previous  attenuation  derived 
data  obtained  during  wind  speed  6  knots. 
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+3 
the  ratio  6c/ (c)  behaved  as  p   .   That  deduction  was  based  on  attenuation 

measurements  performed  at  wind  speeds  of  6  knots.   The  slope  is  plotted  in 
Figure  6  for  comparison  with  the  present  data.   It  appears  that  both  the 
24.4  and  95.6  kHz  data  approximately  follow  that  same  depth  dependence 
except  for  the  contrary  slope  due  to  the  previously  mentioned  24.4  kHz 
resonant  bubbles  nearer  the  surface.   It  is  interesting,  also,  to  compare 
the  magnitudes  of  the  bubble  effect,  as  a  function  of  wind  speed.   This  is 
possible  if  one  assumes  that  the  waters  in  BASS  Strait  are  biologically 
similar  to  those  off  San  Diego,  and  that  the  cause  of  the  sound  speed  dif- 
ference is  bubble  presence  due  to  the  greatly  different  wind  speed.   The  ratio 
of  the  differential  sound  speeds  is  approximately  6  to  1,  which  is  close  to 
the  wind  speed  ratio  of  the  two  experiments  (28  knots  to  6  knots).  Therefore, 
the  proportionality  of  volume  fraction  and  differential  sound  speed  to  wind 
speed,  assumed  in  reference  15,  is  roughly  verified  here.  More  data  of  this 
type  are  needed  to  determine  whether  the  simple  proportionality  is  indeed 
correct  and,  if  so,  the  limit  of  this  proposed  law. 

B.   Average  Speed  of  Sound  as  a  Function  of  Sound  Frequency 

The  dispersive  effect  of  ambient  bubbles  at  sea  has  been  predicted  on 

the  basis  of  excess  attenuation  measurements   ,  and  has  been  measured  directly 

11  14 
in  two  previous  experiments   '   .   Although  an  attempt  to  repeat  the  ambient 

experiment  in  BASS  Strait  was  frustrated  by  the  tape  recorder  failure  in 

BASS  II,  data  of  special  interest  with  respect  to  bubbles  near  structures  was 

obtained  during  BASS  I.   The  previously  discussed  electronic  phase  shift, 

measured  in  BASS  II,  and  smoothed  between  24  and  96  kHz  was  used  to  reduce 

the  data.   A  photograph  of  the  frothy  surface  of  the  water  next  to  the  tower 

docking  area  and  about  5  meters  above  the  measuring  equipment  is  shown  in  Fig. 
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Fig  7.   Photograph  of  the  ocean  surface  at  the  base  of  the  tower  during 
the  data  collection  for  the  graph  of  Fig  8.   The  equipment  is  at 
depth  15  feet,  next  to  the  tower. 
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The  sound  speed  dispersion  is  shown  by  solid  line  in  Figure  8.  The  curve  of 
reference  14  is  shown  by  dashed  line.  The  peaks  and  troughs  of  differential 
speed  are  the  greatest  that  we  have  ever  determined. 

C.  Average  Bubble  Fraction 

The  estimated  asymptotic  differential  speed  of  -15  meters/sec,  of 
Figure  8,  when  used  in  conjunction  with  equation  (11)  of  section  3,  leads  to 
the  value  of  the  average  total  fractional  volume  of  gas  in  bubble  form  to 
water,  U  =  1.2  x  10   for  the  region  near  the  tower.   This  value  is  approxi- 
mately 30  times  the  fraction  for  ambient  bubbles  found  at  3.3  m  depth  in  a 
6  knot  wind  condition,  10  times  the  values  at  depths  to  7.6  m  in  the  28  knot 

wind  (Table  2)  away  from  the  tower  in  BASS  II,  and  6  times  the  value  measured 

4 
directly  by  microscope  near  breaking  waves  at  shore  .   One  might  suspect  that 

the  optical  measurement,  since  it  is  effective  only  for  the  larger  bubbles, 
underestimates  the  total  fraction.   It  is  interesting  to  observe  that  the 
peaks  of  sound  speed  are  at  frequencies  close  to  those  found  in  previous  work 
although  the  experimental  frame  and  particular  hydrophones  and  preamplifiers 
were  different  in  the  different  experiments,  partly  to  eliminate  the  possi- 
bility of  instrumental  artifacts. 

D.  Fluctuations  of  Bubble  Fraction  and  Sound  Speed 

Fluctuations  in  sound  speed  may  be  caused  by  fluctuations  in  temperature, 
salinity,  depth,  fraction  of  gas  in  bubble  form,  resonant  frequencies  of 
bubbles  or  distribution  of  bubbles  with  radius.   From  a  consideration  of 
Wilson's  equation  for  the  speed  of  sound  in  a  bubble-free  region,  it  is  clear 
that  salinity  and  depth  variations  (due  to  change  of  sea  surface  displacement 
with  time)  are  generally  quite  unimportant  compared  to  the  temperature  varia- 
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tions.   It  is  of  interest,  however,  to  compare  the  fluctuations  of  sound 
speed  due  to  temperature  variation  with  those  attributable  to  the  different 
bubble  effects.   This  can  be  done  by  means  of  equation  (19). 

Table  3  displays  the  results  of  calculations  of  fluctuations  during  the 
BASS  II  depth  experiment.   The  experiment  was  conducted  from  1110  to  1255 
(24.4  kHz)  and  1255  to  1442  (95.6  kHz)  on  Saturday,  17  February  1973.   The 
equipment  was  lowered  to  five  depths,  3.0,  4.6,  6.1,  7.6,  13.7  m,  at  a  loca- 
tion about  100  feet  from  the  tower,  which  was  evidently  free  of  tower  effects, 
The  hydrophone  separation  was  0.99  m.   The  wind  speed  of  25-30  knots  measured 
at  height  37  m  above  the  ocean  surface  leads  to  a  predicted  rms  height  of 
waves,  <3TT  =  1.1  m,  assuming  the  Pierson-Moskovitz  spectrum  for  the  fully 
developed  sea.   This  wind  had  persisted  for  about  24  hours  prior  to  the 
acoustics  experiment.   A  swell  wave  of  estimated  peak  to  trough  height  of  2 
meters  appeared  to  be  present.   Under  these  conditions  the  rms  variation  of 
the  average  speeds  of  sound  at  24.4  kHz  and  95.6  kHz  at  five  depths  are  given 
in  column  two  of  Table  3.   Column  four  presents  the  calculated  rms  values  of 
c  which  is  calculated  from  the  variance  of  the  temperature  and  the  constant, 
b,  of  the  linearized  equation  (15).   The  variance  of  the  excess  speed  due  to 
bubbles  is  obtained  by  subtraction  (equation  (19))  and  its  square  root,  the 
rms  value,  is  presented  in  column  five,  "EXPT".   It  is  seen  to  be  greater 
than  the  rms  value  due  to  temperature,  (see  footnote  Table  3)  except  for  95.6 
kHz  at  the  greatest  depth  of  the  experiment,  13.7  m,  where  we  can  assume  that 
the  bubble  population  is  negligible.   Since  variances  add,  it  is  appropriate 

to  compare  the  variance  of  sound  speed  due  to  temperature  fluctuations, 

2 
Var  c  =  G        ,  as  a  percentage  of  the  total  variance  of  sound  speed  fluctua- 

co 
tions,  as  shown  in  column  six.   The  relative  variance  due  to  temperature 
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fluctuations  is  particularly  small  at  the  frequency  24.4  kHz,  and  particu- 
larly large  at  both  frequencies  at  the  greatest  depth  due  to  the  deficiency 
of  bubbles  there.   The  distinction  between  the  two  frequencies  will  be  dis- 
cussed later  after  we  have  looked  at  the  temporal  correlations. 

Finally,  the  standard  deviation  of  fluctuations  of  the  bubble  fraction 
is  shown  in  the  last  column  for  the  24.4  kHz  sound,  where  a  calculation  is 
possibly  by  using  equation  (20).   The  standard  deviations,  0~ '    ,  range  from 
one-third  to  one-eighth  of  the  average  bubble  fractions  themselves,  U,  given 
in  Table  2. 

If  one  has  the  situation:  (a)  the  variance  of  sound  speed  due  to  temper- 
ature fluctuations  is  negligible  and  (b)  the  variance  of  sound  speed  is 
principally  due  to  perturbation  of  the  resonant  frequency  of  bubbles,  it  is 
possible  to  use  equation  (28)  to  estimate  the  rms  value  of  the  excess  speed 
of  sound  fluctuations.   Assumption  (a)  is  verified  by  column  six  of  Table  3, 
except,  of  course,  for  the  greatest  depth,  13.7  meters.   Assumption  (b) 
would  be  applicable  only  if  we  are  considering  a  sound  frequency  at  which 
there  is  a  very  large  population  of  bubbles.   Such  a  situation  can  be  iden- 
tified either  by  a  large  excess  absorption,  a  large  excess  speed  of  sound, 
or  a  large  cross  correlation  between  the  sound  phase  and  the  surface  wave 
height,  such  as  was  observed  in  reference  14.   The  cross  correlation  test 
is  the  best  criterion  to  use  but,  unfortunately,  surface  wave  heights  were 
not  recorded  in  this  experiment.   In  section  F  we  will  show  that  at  95.6  kHz 
the  phase  fluctuation  spectrum  is  highly  responsive  to  the  surface  wave 
height  spectrum.   This  is  pretty  good  evidence  of  a  strong  bubble  population 
resonant  to  95.6  kHz,  so  we  have  proceeded  on  that  assumption. 

The  parenthetical  inserts  of  Column  5  of  Table  3,  presents  the  results 
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of  our  attempt  to  use  the  simple  linear  theory  (equation  (28))  to  predict 
a  k  ,  the  rms  value  of  the  speed  fluctuations  due  to  bubbles  for  the  fre- 
quency 95.6  kHz.   Comparison  with  the  experimental  values  of  column  5  shows 
that  our  predictions  are  an  order  of  magnitude  too  large.   It  is  not  dif- 
ficult to  find  the  reason  for  this  overestimation.   The  flaw  is  that  we 
have  used  our  simple  theory  (equations  21-28)  which  was  developed  for  a 
single  radius  population  with  A  <  0.05  in  attempting  to  predict  frequency 
swings  that  were  as  great  as  A  =  0.16  in  a  multiradius  population  in  the 
high  seas  of  BASS  II.   The  theory  overestimates  the  magnitude  of  the  speed 

fluctuations  because  it  uses  the  slope  near  f  =  f  for  the  calculation, 

R 

when,  in  fact,  this  slope  is  the  maximum  slope  and  is  applicable  only  for 
bubbles  of  precisely  f  =  fR.   For  bubbles  slightly  off  resonance  and  for 
larger  excursions  of  pressure,  the  changes  are  very  much  less.   Further- 
more, for  very  large  excursions  in  pressure,  the  phase  shift  becomes 
"counterproductive"  because  the  sections  of  the  dispersion  curve  over  the 
hump  start  to  produce  phase  shifts  that  cancel  the  effect  near  f  =  f. 
We  believe,  therefore,  that  the  order-of -magnitude  disagreement  between  the 
experimental  and  theoretical  parts  of  column  5  are  not  a  contradiction  of 
our  theory  of  the  probable  cause  of  the  speed  fluctuations  at  95.6  kHz. 
Additional  confirmation  follows  in  section  F. 

E .   Temporal  Correlations 

The  tentative  conclusions  about  the  acoustical  effect  of  bubbles,  which 
we  obtained  by  consideration  of  the  variances,  are  reinforced  when  we  con- 
sider the  temporal  correlations.   Consider  Figure  9,  which  shows  the  tem- 
poral correlations  of  the  temperature  at  the  four  depths,  3.0,  4.6,  6.1 
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and  13.7  m  during  the  24.4  kHz  run.   The  gross  appearance  is  of  an  expo- 
nential or  Gaussian  correlation  function  with  a  periodicity  of  about  6 
seconds  (see  reference  18) .   This  corresponds  to  a  narrow  band  noise 
centered  around  a  peak  at  0.16  Hz,  which  is  approximately  the  frequency  at 
which  peak  energy  density  is  found  for  the  wind  wave  system  that  existed 
at  the  time  of  the  experiment.   It  is  noted  that  the  correlation  period  is 
a  little  shorter  (higher  frequency)  for  the  temperature  fluctuations  at 
3.0  m  depth  compared  to  the  greater  depths;  this  would  be  accounted  for  by 
filtering  action  due  to  depth.   That  is,  at  greater  depths,  the  frequency 
spectrum  of  the  particle  velocities  associated  with  the  surface  wave  system 
will  be  shifted  toward  a  lower  peak  frequency  because  of  the  greater  attenu- 
ation and,  therefore,  the  more  effective  removal  of  high  frequencies.   Super- 
imposed on  the  large  period  variation  is  a  correlation  periodicity  of  about 
1  second.   We  believe  that  this  fine  structure  is  an  artifact  produced  by 
oscillation  of  the  cable  and/or  the  frame  which  held  all  of  the  equipment. 
The  fine  structure  is  repeated  in  all  the  records  of  temperature,  components 
of  velocity,  and  phase  at  13.7  m.  We  believe  that  this  fluctuation  exists 
also  in  the  phase  at  other  depths  but  is  overwhelmed  by  the  bubble  effects 
that  we  were  studying  at  nearer  surface  positions. 

Turning  now  to  the  temporal  autocorrelations  of  the  sound  phase,  we 
consider  Figures  10  and  11,  for  24.4  and  95.6  kHz,  respectively.  In  Figure 
10,  we  see  that,  at  the  greatest  depth,  13.7  m,  the  temporal  correlation  of 
the  sound  phase  at  24.4  kHz  is  almost  a  replica  of  the  temperature  correla- 
tions of  Figure  9.  The  fine  structure  appears  with  the  same  periodicity  of 
about  1  second,  and  the  6  second  periodicity  due  to  motion  of  the  surface 
wave  system  is  clearly  seen  as  well.   Recalling  our  analysis  of  the 
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variances,  this  is  a  case  in  which  the  temperature-caused  variance  of  the 
speed  is  greater  than  the  bubble-caused  variance.   It  is,  therefore,  not 
surprising  to  find  the  sound  phase  correlation  resembling  the  temperature 
correlation  at  24.4  kHz  as  the  temperature  microstructure  moves  with  respect 
to  the  frame.   There  is  only  a  hint  of  the  fine  structure  at  the  depth  7.6 
meters  where  Table  3  has  shown  us  that  the  variance  of  sound  speed  fluctua- 
tions due  to  temperature  is  only  0.2%  of  the  variance  due  to  bubble  effects. 
At  shallower  depths  there  is  even  less  fine  structure  because,  even  though 
the  frame  moved,  the  phase  fluctuations  were  due  to  bubbles  rather  than  to 
temperature  blobs.   It  is  our  belief  that  the  decorrelation  of  the  sound 
phase  at  all  depths  except  13.7  m  is  principally  a  measure  of  the  random- 
ness of  the  strong  bubble  populations  at  these  nearer-surface  positions. 
There  is  virtually  no  response  to  the  ocean  wave  system  at  the  near-surface 
positions  in  the  case  of  the  24.4  kHz  sound. 

Figure  11,  the  temporal  autocorrelation  of  sound  phase  fluctuations  at 
frequency  95.6  kHz  presents  a  completely  different  story.   Only  at  13.7  m 
is  there  a  hint  of  the  temperature -mot ion  fine  structure.   The  predominant 
general  feature  is  the  narrow  band  noise  character  of  the  correlation  func- 
tion with  a  periodicity  of  approximately  6  seconds,  that  is,  the  frequency 
0.16  Hz  again,  the  center  of  the  surface  wave  energy  spectrum.  We  propose 

that  the  explanation  of  these  correlation  functions,  which  are  so  close  to 

23 
what  is  to  be  expected  of  a  surface  wave  correlation  function   ,  is  due  to 

the  responsiveness  of  the  resonant  frequency  of  predominant  bubble  popula- 
tions to  surface  wave  pressure  variations  as  explained  in  Section  III,  B, 
and  as  used  in  Section  VI,  D,  to  calculate  the  theoretical  variance  in 
Column  5  of  Table  3.  We  do  not  have  the  directly  measured  bubble  distribu- 
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tions  to  prove  absolutely  that  there  was  a  large  bubble  population  resonant 
at  95.6  kHz  at  all  depths.   However,  it  should  be  recalled  that  the  resonant 
frequency  is  proportional  to  the  5/6  power  of  the  ambient  pressure.   This 
means  that  a  large  bubble  population  resonant  at  frequency  f   at  depth  3.0  m 
will  be  resonant  at  2fQ  at  18  m.   Theref6re,  it  should  not  be  surprising  to 
find  effects  of  large  bubble  populations,  resonant  at  95.6  kHz,  at  all  depths 
of  this  experiment  conducted  in  a  sea  churned  up  by  25-30  knot  winds.   In 
fact,  there  is  ample  ancillary  evident  of  a  large  population  of  bubbles  at 
95.6  kHz,  as  will  be  seen  in  Section  F. 

Two  temporal  cross  correlations  have  been  plotted  in  an  attempt  to 
further  characterize  the  effect  of  the  medium  on  the  sound  phase.   Figure  12 
shows  the  cross  correlation  of  flow  velocity  fluctuations  (the  axial  compo- 
nent of  the  velocity  of  the  water  with  respect  to  the  frame  and  sound  field 
axis)  and  24.4  kHz  sound  phase  variations.   It  is  observed  that  the  long 
term  periodicity  of  6  seconds  is  again  found  for  both  3.0  m  and  13.6  m 
depths,  but  the  fine  scale  fluctuation  of  approximately  1  second  periodicity 
is  clearly  identifiable  only  at  the  greater  depth.   The  fact  that  the  fine 
structure  is  almost  non-existent  on  the  correlogram  for  3.0  m  is  presumably 
due  to  the  very  much  larger  peak  frequency  orbital  flow  velocity  components 
which  will  tend  to  swamp  the  other  temperature-caused  fluctuations  of  sound 
phase  variation  after  normalizing.   Figure  13  is  a  correlogram  of  the 
temporal  correlation  of  flow  with  temperature.   It  is  almost  identical  to 
the  previous  figure  for  the  same  reasons. 

It  is  observed  that  the  correlation  envelopes  at  95.6  kHz  require  much 
more  time  to  decay  than  their  counterparts  at  24.4  kHz,  particularly  the 
shallow  depth  correlograms  at  24.4  kHz.   The  interpretation  that  we  make  is 
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that  there  is  greater  variability  in  the  instantaneous  total  number  of 
bubbles  relative  to  their  average  number  than  there  is  variability  in 
temperature.   Tables  2  and  3  show  that  CTT,/U  ranges  from  about  one-third 
to  one- tenth  as  depth  increases.   Notice  that  the  near-surface  24.4  kHz 
correlograms  are  down  to  1/e  of  their  value  in  about  1  second,  whereas 
the  envelope  of  the  temperature-dominated  24.4  kHz  correlogram  at  13.7  m 
depth  is  still  above  1/e  after  as  much  as  7  seconds'  delay.   If  this 
reasoning  is  correct,  the  spatial  correlation  function  of  number  of 
bubbles  (if  this  could  be  measured)  would  have  a  smaller  correlation  length 
than  the  spatial  correlation  function  of  excess  temperature. 

Using  the  same  reasoning  as  above,  we  must  conclude  that  the  number  of 
resonant  bubbles  at  95.6  kHz,  which  determines  the  character  of  those  cor- 
relograms, must  be  reasonably  constant  over  a  7  second  lag  time  (the  maxi- 
mum calculated).   Since  the  95.6  kHz  bubbles  are  the  smallest  ones  measured 

in  this  experiment,  it  appears  that  the  small  bubbles  retain  their  identity, 

Aa 

—  ^  ^r/2  «*  .05,  and  their  gross  locations,  for  a  time  of  at  least  7  seconds. 

F.   Spectra 

it 

Data  processing  of  the  spectra  of  0,  T  and  v  have  been  done  with  dif- 
ferent time  lags  and  different  resolutions  and  band  width  in  order  to  clarify 
the  spectral  implications  already  deduced  from  the  correlation  functions. 
The  previously  described  "fine  structure"  that  is  so  prominent  in  the  temper- 
ature correlations  (Figure  9)  and  elsewhere,  turns  out  to  be  an  harmonic 
series  with  0.5  Hz  fundamental  and  strong,  even,  harmonics  particularly  at 

*To  avoid  confusion  with  the  symbol,  u,  used  for  volume  fraction,  we  will 
use  v  for  the  horizontal  material  velocity  in  the  direction  of  the  sound 
beam. 
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1,  2,  and  4  Hz.   These  components  are  strongest  in  the  material  velocity  and 
temperature  spectra,  and  we  believe  that  they  represent  motion  of  the  equip- 
ment frame.   Since  this  oscillating  motion  is  in  the  presence  of  much  larger 
D.C.  velocities  (e.g.,  v  -  0.64  m/sec  and  o     -   0.064  m/sec  in  one  case  -  see 
Appendix  3)  the  even  harmonics  are  emphasized.   Further,  since  the  oscillat- 
ing motion  of  the  frame  takes  the  temperature  sensors  to  positions  of  dif- 
ferent temperature  (due  to  the  vertical  temperature  gradient),  the  temper- 
atures show  similar  spectra.   In  the  cases  where  temperature  effects,  rather 
than  bubble  effects,  dominate  the  fine  structure  is  evident  also  in  the 
sound  phase  fluctuations.   Since  the  fine  structure  is  an  artifact  that  is 
characteristic  of  our  mounting  structure  rather  than  the  ocean,  we  will 
not  discuss  this  aspect  of  the  spectra  any  further. 

A  higher  resolution  (F*  =  0.025  Hz),  smaller  bandwidth  (Fmax  =  0.6  Hz), 

presentation  of  the  horizontal  velocity  component  of  material  velocity  at 

19 
five  depths  is  shown  in  Figure  14.   Previous  study  '  has  shown  that  the 

spectral  components  of  the  material  velocity  in  the  ocean  are  well  corre- 
lated with,  and  representative  of,  the  ocean  surface  height  spectrum, 
attenuated  and  filtered  by  the  depth  dependent  factor  exp(-Kz),  where  K  is 
the  wave  number  of  the  spectral  component  of  the  orbital  velocity.   The 
velocity  spectra  of  Figure  14,  do,  indeed,  resemble  those  of  an  ocean  sur- 
face height  spectrum  (confirmed  also  in  Figure  12  of  the  BASS  II  Trials 
Officer's  Report4^)  with  peak  power  spectral  density  (PSD)  at  the  frequency 
0.175  Hz. 

The  behavior  of  the  PSD  of  $V(F)  at  the  peak  frequency  (0.175  Hz)  has 

*To  avoid  confusion,  we  will  use  the  upper  case  "F"  for  modulation 
frequency  and  the  lower  case  "f"  for  sound  carrier  frequency. 
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Fig  14.   Power  spectral  densities  of  fluctuations  in  the 
horizontal  material  velocity  at  3.0,  4.6,  6.1, 
7.6  and  13 . 7  m  during  the  95.6  kHz  sound  phase 
study.   BASS  II. 
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been  identified  and  replotted  in  Figure  15  as  a  function  of  depth  in  order 
to  compare  with  theory.   At  this  frequency,  the  surface  wave  number  is 
K  =  0.116  m   .   If  the  spectral  component  is  following  the  attenuation  law 
with  depth,  20  log-^Q  of  exp(-Kz)  gives  a  change  of  level,  with  depth  of 
1.01  dB/m.   A  best  line  of  this  slope  has  been  drawn  through  the  five  data 
points  available,  and  it  is  seen  to  quite  well  represent  the  behavior  of  the 
observed  material  velocity  component.  We  conclude  that  the  peak  component 
of  the  underwater  particle  velocity  spectrum  in  BASS  II  was  principally  due 
to  the  surface  wave  height  spectrum.   The  point  at  which  the  higher  frequency 
components  of  the  velocity  spectrum  are  dominated  by  turbulent  (non-orbital) 
components  is  not  clear. 

The  spectra  of  the  temperature  fluctuations  at  five  depths  during  the 
95.6  kHz  run  are  plotted  in  Figure  16.   Since  the  temperature  spectrum  in- 
volves the  manner  in  which  a  spatially  distributed  temperature  microstruc- 
ture  is  convected  past  the  sensor  by  the  complex  motion  of  the  medium  with 
respect  to  the  sensor,  the  spectrum  is  much  more  difficult  to  comprehend 
than  the  velocity  spectrum.   Nevertheless,  some  understanding  can  be  deduced. 
For  example,  the  spectra  are  due  to  rates  of  change  of  the  temperature  and 
are,  therefore,  measured  by  the  material  change, 

The  local  rate  of  change,  d/dt,  can  be  assumed  to  be  negligible  com- 
pared to  the  convective  change  if  the  Taylor  Hypothesis  of  frozen  convected 
scalars  is  accepted.   Since  vertical  temperature  gradients  are  more  likely 
than  horizontal  ones,  the  convective  change  will  be  maximum  in  the  z  direc- 
tion, that  is,  the  term  w  BT/dz.   Unfortunately,  only  one  electromagnetic 
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material  velocity  fluctuations  at  the  peak  frequency 
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flowmeter  was  available  and  it  was  set  to  measure  the  horizontal  components 
rather  than  the  vertical.   However,  our  previous  discussion  of  Figure  14 
suggests  that  the  component  at  0.17  Hz  was  at  the  peak  of  both  the  pressure 
and  the  velocity  spectrum,  and,  therefore,  that  the  z  component  of  the 
orbital  velocity,  w,  would  also  peak  at  this  frequency.   Now  consider  Fig- 
ure 16.  We  see  that  even  though  the  vertical  component  of  the  material 
velocity  is  greater  at  3.0  than  at  4.6  m,  the  temperature  fluctuation  spec- 
trum at  the  peak  frequency  is  less  at  the  shallower  position.   This  can  only 
be  due  to  a  vertical  temperature  microstructure  that  has  a  large  gradient 
BT/Bz  at  the  lower  level.   This  reversal  of  the  expected  decrease  in  $m  with 
increasing  depth  does  not  occur  at  any  other  level,  at  the  peak  frequency. 
An  inversion  at  approximately  5  m  in  the  BT  at  1010  and  a  "twitch"  at  about 
4  m  in  the  BT  at  1415  (see  Appendix  2)  support  our  observation  that  the 
microstructure  at  4.6  m  was  distinctive.   The  effect  also  shows  in  the  C^ 
column  of  Table  3,  during  the  95.6  kHz  run. 

Another  feature  of  the  temperature  spectra  is  that  the  prominent  peaks 
of  energy  at  0.17  and  0.27  Hz  at  the  shallower  depths  appear  to  be  shifting 
to  higher  frequencies  as  the  depth  is  increased.   This  effect  appears  also 
in  the  velocity  spectra,  and  may  be  evidence  of  underwater  turbulence,  uncor- 
rected with  the  surface  height,  which  is  becoming  increasingly  significant 
as  the  depth  increases.   If  this  is  the  true  explanation,  it  appears  that 
the  turbulence  spectrum  has  energy  peaks  at  higher  frequencies  than  those 
of  the  orbital  velocity  spectrum.   There  is  some  confirmation  of  this  specu- 
lation in  the  greater  than  expected  value  of  the  velocity  at  13.6  m  in 
Figure  15.   The  temperature  spectrum  appears  to  imitate  these  peaks  in  the 
velocity  spectrum  and  therefore  to  be  behaving  somewhat  like  the  passive 
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20 
convected  scalar  of  Batchelor's  theory 

Spectra  of  the  fluctuations  of  sound  phase  are  shown  in  Figures  17,  18, 
and  19.   Compared  to  95.6  kHz,  the  spectra  for  sound  frequency  24.4  kHz  are 
generally  flatter  and  closer  to  a  Gaussian  form,  as  would  be  expected  from 
the  near-Gaussian  correlation  functions  (figures  17,  19).   The  broader  spectra 
for  24.4  kHz  are  undoubtedly  a  consequence  of  the  strong  dependence  of  the 
sound  speed  on  the  bubble  fraction,  U,  which  is  likely  to  be  Gaussian  dis- 
tributed (statistically)  and  to  have  a  Gaussian  spatial  correlation  function. 

We  have  already  pointed  out  the  broad  PDF  implied  by  0.12  ^  _H  =s  0.30.   There 

U 

is  essentially  no  evidence  of  the  peak  surface  wave  frequency  (0.17  Hz)  in 
the  24.4  kHz  spectra. 

The  predominant  feature  of  the  phase  fluctuations  for  95.6  kHz  is  the 
strong  peak  at  the  ocean  surface  wave  peak  frequency  of  0.17  Hz*.   Figures 
18  and  19  show  this  peak  to  be  common  to  all  depths  and  to  generally  decrease 
in  magnitude  with  increasing  depth;  this  provides  additional  evidence  of  our 
explanation  that  the  fluctuations  at  this  frequency  are  predominantly  due  to 
change  of  radius  caused  by  the  surface  wave  pressure.   Away  from  this  domi- 
nant frequency,  the  spectra  at  depths  3.0,  6.1,  7.6  and  13.7  m  are  almost 
indistinguishable.   Beyond  F  =  0.3  Hz,  the  95.6  kHz  spectrum  at  depth  4.6  m 
shows  larger  PSDs  than  at  3.0  m  (as  does  the  phase  spectrum  for  frequency 
24.4  kHz);  this  is  apparently  due  to  the  larger  temperature  fluctuations  at 
4.6  m  previously  observed  in  Figure  16.   Except  for  this  anomaly  at  4.6  m, 


*For  these  graphs  the  resolution  is  0.05  Hz  so  that  what  was  formerly 
identified  as  a  peak  at  0.175  Hz,  now  is  found  at  0.15  rather  than  0.20  Hz; 
in  all  cases  we  have  compromised  and  arbitrarily  assigned  the  unresolvable 
value  0.17  Hz  to  this  peak  frequency  so  that  we  will  not  cause  confusion  in 
the  identification  by  the  reader. 
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for  F  >  0.17  Hz,  the  power  spectral  density  of  the  phase  fluctuations  for 

95.6  kHz  sound  follows  precisely  the  F   law  that  would  be  expected  for  a 

21  22 
surface  wave  pressure  effect   '  '  for  a  fully-developed  sea.   This  is  shown 

most  dramatically  in  Figure  20,  where  we  have  replotted  some  of  the  data  of 
Figures  18  and  19  on  log-log  paper  in  order  to  show  the  power  law  dependence. 
At  low  frequencies  the  phase  fluctuation  spectrum  is  almost  a  twin  to  the 
ocean  surface  wave  spectrum  and  the  fifth  power  law  for  F  >  F_  is  found  at 
all  depths  for  a  range  of  frequencies  from  just  past  F  =  0.17  Hz  to  approxi- 
mately F  =  0.5  or  0.7  Hz.   The  exception,  again,  is  the  case  for  depth  4.6  m 
where  the  fluctuating  temperature  effect  begins  to  take  over  from  the  pres- 
sure dependence  for  F  >  0.35  Hz  (graphs  for  4.6  m,  6.1  m,  and  13.7  m  not 
shown  to  avoid  crowding).   The  displaced  lines  of  slope  F   have  been  re- 
plotted  versus  depth  (not  shown).   The  variation  of  the  spectral  level  with 
depth  has  been  determined  to  be  approximately  exp(-H/4.7  m) ,  which  compares 
fairly  well  with  a  determination  of  (n/a)da  from  decay  of  attenuation  with 
depth  (reference  10),  that  behaved  as  exp(-H/6.7  m) . 

As  a  conclusive  check  to  show  that  the  spectrum  of  phase  fluctuations 
at  95.6  kHz  was  not  determined  by  the  particle  velocity,  the  phase  shift  due 
to  spectrum  of  underwater  material  velocity  components,  $(F)  =  Ae       F   , 
was  plotted  separately  and  compared  with  the  data  for  the  3.0  m  spectrum  of 
phase  fluctuations  (not  shown);  the  velocity  theory  did  not  fit  the  data. 
We  conclude  that  the  phase  modulation  at  95.6  kHz  is  due  to  changes  of  radius 
of  resonant  bubbles  in  response  to  pressures  caused  by  the  major  frequency 
components  of  the  surface  wave  height  spectrum. 
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VIII.   APPENDICES 

APPENDIX  1 

BASS  I  EXPERIMENTS 

I.   Data 

1.  Fluctuations  experiment  at  separations  1.047,  2.045,  4.046  m 

2.  Decorrelation  rate:  data  at  15,  25,  35  ft  depth.   One  separa- 

tion. 62.15  kHz. 

3.  Bubble  sweep:   data  at  28  frequencies  from  15  to  97  kHz,  15  ft 

depth. 

BASE  II  EXPERIMENTS 

I.   Data 

A.   Acoustical  data 

1.  Frequency  run:  four  minute  runs  from  10  kHz;  single  wavelength 
differences  from  50  to  75  kHz,  selected  frequencies  elsewhere. 
Depths:   10,  70  (max),  20,  10  ft.   Estimated  Duration  of  Record- 
ing (EDR) :  4  hours;  Estimated  Duration  of  Experiment  (EDE)  : 

6  hours. 

2.  Depth  run:  10  minute  runs  at  two  frequencies,  at  six  depths. 

a.  Set  frequency  at  25  kHz*, vary  depth  10,  15,  20,  25,  50, 
and  70  feet. 

b.  Set  frequency  at  95  kHz,  vary  depth  70,  50,  25,  20,  15, 
and  10  feet;  EDR:  2  hours;  EDE:  3  hours 

3.  Diurnal  run:  sound  speed  as  function  of  time 


a.   Set  depth  10  feet  below  trough.  Consecutive  10  min  runs  at 
2  5  kHz   and  65  kHz;  repeat  every  hour  for  24  hours. 
EDR:  8  hours;  EDE  24  hours 

^Frequency  25  kHz  is  frequency  for  which  S /N  is  at  least  35  dB 
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4.   Wind  run:  sound  speed  during  changing  wind  conditions 


a.   Set  depth  10  feet  below  trough.   Do  continual  recording  at 
25  kHz   during  rapidly  chan; 
EDR:  3  hours;  EDE:  3  hours. 


25  kHz   during  rapidly  changing  wind  conditions. 


•k 

5.   Range  run;  sound  speed  as  a  function  of  range;  frequencies  25 
and  95  kHz  at  four  ranges  (10  min  runs)  at  depth  10  feet, 
then  repeat  at  70  feet. 
EDR:  3  hours;  EDE:  4  hours 

B.  Auxiliary  data  needed  during  all  NPS  experiments 

1.   Two  horizontal  components  of  flow  velocity  (u  and  v) ;  temper- 
ature, T;  velocimeter  fluctuations  and  frequency  averages. 

C.  Hourly  observations  needed  during  NPS  experiments 

1.  XBT  from  tower 

2.  Meteorological  observations:   wind  direction,  speed,  sky 
coverage,  precipitation  (height  of  anemometer) 

3.  Oceanographic  observations:  sea  state  and  swell  height,  and 
periodicity  and  direction,  white  caps,  foam. 

II.   Data  Processing 

A.  Time  series  graphs  to  determine  validity  of  recordings 

B.  Mean  values  of  all  data 

C.  Variance  wrt  mean  value  (=  zero)  for  all  data 

D.  Temporal  Correlations:   corrected  c,  all  data  of  depth  run  (6  depths, 
2  frequencies)  and  six-10  min  extracts  of  wind  run  v,  T:  one  depth 
run  and  six-10  min  extracts  of  wind 

JL. 

E.  Auto  spectra:  corrected  c  at  10  and  70  feet;  at  25  and  95  kHz  plus 
6  winds,  u,v,T:  at  10  and  70  feet  plus  u  ,  only  for  6  winds 

F.  Cross  Spectra,  Coherences  and  Phase  Angles:  for  same  data  as  in  F . 

NOTES:   Speed  of  sound,  c,  referred  to  above,  is  always  speed  of  sound 
corrected  for  phase  shift  due  to  longitudinal  velocity,  u. 

60. 


APPENDIX    2 


BATHY  THERMOGRAMS 
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APPENDIX  3 


MATERIAL  VELOCITIES  CALCULATED  FROM  ELECTROMAGNETIC 
FLOWMETER  OUTPUT  DURING  BASS  II  -  17  FEB  1973 

u  =  horizontal  speed  in  direction  of  sound  path 
v  =  horizontal  speed  perpendicular  to  sound  path 


Clock 
Time 

Depth 
m 

Concurrent 
sound  freq. 
kHz 

<u> 
m/sec 

m/sec 

<v> 
m/sec 

m/sec 

1110 

3.0 

24.4 

0.23 

0.15 

0.33 

0.18 

1235 

13.7 

24.4 

0.64 

0.064 

- 

— 

1258 

13.7 

95.6 

0.64 

0.071 

- 

- 

1318 

7.6 

95.6 

0.61 

0.08 

- 

- 

1343 

6.1 

95.6 

0.56 

0.09 

- 

- 

1401 

4.6 

95.6 

0.51 

0.10 

- 

- 

1430 

3.0 

95.6 

0.53 

0.11 

— 

— 

NOTE:   v  component  output  failed  after  1130. 
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